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The Main Injector Particle Physics Experiment (MIPP FNAL E-907) at 
Fermilab - status and plans 

Rajendran Raja a 

a Fermilab, P.O. Box 500, Batavia, ILlinois 60510, U.S.A 



We describe the status of the Main Injector particle production Experiment (MIPP) at Fermilab which has 
to date acquired 18 million events of particle interactions using (5 GeV/c-120 GeV/c) tt ± ,K ± and p ± beams on 
various targets. We describe plans to upgrade the data acquisition speed of MIPP to make it run 100 times faster 
which will enable us to obtain particle production data of unprecdented quality and statistics on a wide variety 
of nuclear targets including nitrogen which is of importance to cosmic ray physics. 



1. Current Status of the MIPP Experi- 
ment 

We describe status report on the MIPP pQ ex- 
periment and its performance to date. MIPP is 
situated in the Meson Center beamline at Fermi- 
lab. It received approval [2 J in November 2001 
and has installed and operated both the experi- 
ment and a newly designed secondary beamline in 
the interim. It received its first beams in March 
2004, had an engineering run to commission the 
detector in 2004 and had its physics data-taking 
run in the period January 2005-March 2006. The 
experiment is currently busy analyzing its data. 

MIPP is designed primarily as an experiment 
to measure and study in detail the dynamics asso- 
ciated with non-pert urbative strong interactions. 
It has nearly 100% acceptance for charged parti- 
cles and excellent momentum resolution. Using 
particle identification techniques that encompass 
dE/dx, time-of-flight [3], Multi-Cell Cerenkov [4] 
and a Ring Imaging Cerenkov (RICH) detector [5] , 
MIPP is designed to identify charged particles at 
the 3a or better level in nearly all of its final state 
phase space. MIPP has acquired data of unpar- 
alleled quality and statistics for beam momenta 
ranging from 5 GeV/c to 90 GeV/c for 6 beam 
species (tt ± 1 K ± and p ± ) on a variety of targets 
as shown in Figure [TJ 

An important aspect of MIPP data-taking was 
the measurement of particle production off the 
NuMI [6J target in order to minimize the sys- 
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Acquired Data by Target and Beam Energy 
Number of events, x 10 6 
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Figure 1. The data taken during the first MIPP 
run as a function of nucleus. The numbers are in 
millions of events. During the last month of the 
run, the Jolly Green Giant magnet coils devel- 
oped shorts. This time was used to acquire data 
without the TPC for exploring the feasibility of 
measuring the charged kaon mass using the RICH 
radii. 



tematics in the near /far detector ratio in the MI- 
NOS [6 J experiment. MIPP also made measure- 
ments with proton beams off various nuclei for 
the needs of proton radiography [2J. 

Another physics motivation behind MIPP is to 
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restart the study of non-perturbative QCD in- 
teractions. Currently available data are of poor 
quality, and sparsely populate the beam momen- 
tum, pt : and atomic weight phase space that 
makes comparisons between different experiments 
difficult. The MIPP TPC [7J digitizes the charged 
tracks in three dimensions, obviating the need for 
track matching across stereo views. Coupled with 
the particle identification capability of MIPP, the 
data from MIPP would add significantly to our 
knowledge base of non-perturbative QCD. This 
would help test inclusive scaling relations and also 
scaling nuclear reactions. 

1.1. Experimental Setup 

We designed a secondary beam [8] specific to 
our needs. The resonantly extracted protons 
from the Fermilab Main Injector are transported 
down the Meson Center line. They impinge on 
a 20 cm long copper target producing secondary 
beam particles. This target is imaged onto an 
adjustable momentum selection collimator which 
controls the momentum spread of the beam. This 
collimator is re-imaged on to our interaction tar- 
get placed next to the TPC. The beam is tracked 
using three beam chambers and identified using 
two differential Cerenkovs [9 filled with gas, the 
composition and the pressure of which can be var- 
ied within limits depending on the beam momen- 
tum and charge. 

Figure [2] shows the layout of the apparatus. 
The TPC sits in a wide aperture magnet (the 
Jolly Green Giant) which has a peak field of 0.7 
tesla. Downstream of the TPC are a 96 mir- 
ror multi-cell Cerenkov detector filled with C4F10 
gas, and a time of flight system. This is followed 
by a large aperture magnet (Rosie) which runs 
in opposite polarity (at -0.6 tesla) to the Jolly 
Green Giant to bend the particles back into the 
Ring Imaging Cerenkov counter. The RICH has 
CO2 as the radiator and an array of phototubes 
of 32 rows and 89 columns [10]. Downstream of 
the RICH we have an electromagnetic calorime- 
ter [TT]and a hadron calorimeter [12] to measure 
forward-going photons and neutrons. The elec- 
tromagnetic calorimeter provides a means of dis- 
tinguishing forward neutrons from photons and 
will also serve as a device to measure the electron 



content of our beam at lower energies, which will 
be useful for measuring cross sections. 
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Main Injector Particle Production Experiment (FNAL-E907) 




Figure 2. The experimental setup. The picture is 
a rendition in Geant3, which is used to simulate 
the detector. 



MIPP uses dE/dx in the TPC to separate pi- 
ons, kaons and protons for momenta less than « 
1 GeV/c and the time of flight array of coun- 
ters to do the particle identification for momenta 
less than 2 GeV/c. The multi-cell Cerenkov de- 
tector [4] contributes to particle identification in 
the momentum range ^2.5 GeV/c-14 GeV/c and 
the RICH [5 for momenta higher than this. By 
combining information from all counters, we get 
the expected particle identification separation for 
K/p and tt/K as shown in Figure [31 It can be 
seen that excellent separation at the 3<r or higher 
level exists for both K/p and ir/p over almost 
all of phase space. Tracking of the beam parti- 
cles and secondary beam particles is accomplished 
by a set of drift chambers [13] and proportional 
chambers [14] each of which have 4 stereo layers. 

1.2. Some results from Acquired data 

Figure 2] shows the pictures of reconstructed 
tracks in the TPC obtained during the data- 
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taking run. The tracks are digitized and fitted as 
helices in three dimensions. Extrapolating three 
dimensional tracks to the other chambers makes 
the pattern recognition particularly easy. 




mux 




from two different events. 



Figure 3. Particle identification plots for 
pion/kaon separation and for kaon/proton sepa- 
ration as a function of the longitudinal and trans- 
verse momentum of the outgoing final state parti- 
cle. Black indicates separation at the 3<r level or 
better and grey indicates separation at the 1 — 3a 
level. The boxes at largest values of the longitu- 
dinal momenta suffer from lack of kaon statistics. 



Figure [5] shows the distribution of dE/dx of 
tracks measured in the TPC as a function of the 
track momentum in a preliminary analysis of p- 
Carbon data. The TPC is capable of separat- 
ing pions, protons and kaons in the momentum 
range below « 1 GeV/c. Figure [6] shows events 
with multiple rings in the RICH counter. Figure [3 
shows the histogram of ring radii for a +40 GeV 
secondary beam. There is clean separation be- 
tween pions, kaons and protons and their rela- 
tive abundances [15J match expectations. Ap- 
plying the particle identification trigger from the 
beam Cerenkovs enables us to separate the three 
particle species cleanly. The kaons which form 
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Figure 5. Preliminary dE/dx distributions in 
the TPC The scatter plot shows the electron, 
pion, kaon and proton peaks in the distribution 
as a function of the lab momentum for p-Carbon 
data. The second plot is the projection on the 
dE/dx axis for a momentum slice 0.3 GeV/c to 
0.4 GeV/c. 



4% of the beam are cleanly picked out by the 
beam Cerenkov with very simple selection crite- 
ria. These can be made much more stringent with 
offline cuts to produce a very clean kaon beam. 

The ring radius of the particle contains infor- 
mation on the mass of the particle. The pion and 
proton masses are very well known. The charged 
kaon mass, however, currently has measurement 
uncertainties of the order of 60 keV. Improving 
the precision of both charged kaon masses will 
pay dividends in rare K decay experiments in- 
volving charged kaons where the matrix elements 
depend on the kaon mass raised to large pow- 
ers. Towards the end of our physics run, when 
the Jolly Green Giant magnet coils failed, we 
switched off the TPC and acquired data at the 




Figure 6. Examples of events with rings in the 
RICH counter for a 40 GeV/c beam. The x and 
y axes are in cm. 



rate of 300 Hz to investigate how well we can 
measure the charged kaon mass. These events, 
whose statistics are indicated in Figured! are cur- 
rently being analyzed to evaluate the systematics 
involved in such a measurement. 

1.2.1. NuMI target measurements 

MIPP took 1.75 million events using 120 GeV/c 
primary beam protons impinging on the NuMI 
(spare) target. These events will play a cru- 
cial role in the prediction of neutrino fluxes in 
the NuMI beamline and will enable the MI- 
NOS experiment to control the systematics in the 
near /far detector ratios as well as helping them 
understand the near detector performance. Fig- 
ure [8] shows a radiograph of the MIPP measure- 
ments of the MINOS target. The graphite slabs 
and cooling tubes can be seen. These events 
were obtained during the commissioning phase of 
this target measurement where the beam was not 
yet fully focused and aligned on the target. The 
1.75 Million events on the NuMI target were ob- 
tained after the beam was aligned and centered 
on the target. Figure [9] shows the rich ring radii 
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Figure 7. An example of a 40 GeV/c primary 
beam (non-interacting) trigger. The RICH iden- 
tifies protons, kaons and pions by the ring radii. 
The beam Cerenkov detectors can be used to do 
the same. When the beam Cerenkov identifica- 
tion is used, one gets a very clean separation of 
pions, kaons and protons in the RICH. 



NuMI Target Radiograph 




Figure 8. Radiograph of the MINOS target. The 
beam direction is perpendicular to the paper. 



vs momentum of positive tracks originating from 
the NuMI target. Superimposed are the curves 
for known particles. This shows the excellent par- 
ticle identification of the MIPP detector for for- 
ward going particles. 



Figure 9. Preliminary data of RICH ring radii of 
positive tracks from the NuMI target vs momen- 
tum. Superimposed are the expected curves for 
e, /i, K and p particles. 



1.3. Particle production on Nitrogen and 
the question of Cosmic Ray Showers 
in the Atmosphere 

With the upgraded MIPP, we propose to mea- 
sure particle production on a cryogenic nitrogen 
target using positive and negative beams, which 
is needed by experiments measuring cosmic ray 
air showers (Pierre Auger, HiRes etc) and also 
atmospheric neutrinos (Amanda, Ice Cube, Hy- 
perK etc). The problem is illustrated in a recent 
paper [16] which simulates the air showers pro- 
duced by protons of 10 6 GeV energy in the atmo- 
sphere. The shower goes through several genera- 
tions of interactions and produce pions and kaons 
that decay to produce muons and neutrinos. The 
muons and neutrinos are observed in the detec- 
tors and are termed the daughter particles. The 
mesons that produced the muons and neutrinos 
are termed the mother particles and the parti- 
cles that interacted in the atmosphere to produce 
the mother particles are termed the grandmother 
particles in the jargon. 

Figure [10] shows the energy spectrum of the 
grandmother particles (7r, K and p in an air 
shower that are produced by a primary proton 
of 10 6 GeV. The spectrum for the pions peaks at 
100 GeV and the kaons and protons at somewhat 
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higher energies. These particles interact with the 
nitrogen (and oxygen) in the atmosphere to pro- 
duce the atmospheric neutrinos and muons. In 
other words, the beam energies available at MIPP 
are relevant to the simulation of the cosmic ray air 
showers. The muon flux measurement is a criti- 
cal component of estimating the energy scale of 
the cosmic ray shower. MIPP measurements thus 
will help reduce the systematics in the cosmic ray 
energy scale measurements. As the primary cos- 
mic ray energy increases, the peaks in this plot 
do not shift to higher energies. Understanding 
the shower systematics at the peak of this spec- 
trum (i.e MIPP energies) will help the energy sys- 
tematics of cosmic rays of all energies. Figure fTTI 




simulation code QGSJET 01 is used. The sharp 
break in the spectra at 100 GeV is where the two 
codes meet and disagree at places by a factor of 
two. This illustrates the problem. These codes at 
present are "tuned" on single arm spectrometer 
data and disagree with each other. 




E kin (GeV) 



E kin ( GeV ) 



Figure 11. The energy distribution of the grand- 
mother particles in a vertical air shower produced 
by 10 6 GeV proton interacting in the atmosphere 
as a function of distance from the shower center. 
The spectrum peaks at energies relevant to the 
MIPP energy scale. 



Figure 10. The energy distribution of the grand- 
mother particles in a vertical air shower produced 
by 10 6 GeV proton interacting in the atmosphere 
as a function of particle type. These particles in- 
teract further in the atmosphere to produce more 
particles which then decay into muons and neu- 
trinos. The muons are detected at a distance of 
0-500 meters from the shower center at ground 
level. It can be seen that these spectra peak at 
energies relevant to the MIPP energy scale. 



shows the distribution of grandmother particles 
at different lateral distances from the shower cen- 
ter for all particle types. For the lower energy in- 
teractions, the simulation code Gheisha is used to 
simulate the interactions of the particles with the 
atmosphere. For higher energy interactions, the 



Most of the data acquired to date of relevance 
to air showers are over 25 years old and were 
obtained using Beryllium targets and single arm 
spectrometers resulting in a discrete angular cov- 
erage. MIPP will measure the outgoing pion and 
kaon spectrum for proton and pion beams in its 
full secondary beam momentum range. These 
cross sections are the most important in simu- 
lating cosmic ray showers in the atmosphere. In 
addition, it will also have kaon and antiproton in- 
teractions on nitrogen of which virtually nothing 
is known. Being an open geometry experiment, 
the MIPP angular coverage will be continuous, 
not discrete. 

The need for MIPP data is recognized by the 
cosmic ray community, some of whom have joined 
the experiment. 
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1.3.1. MIPP Measurement of tt/K ratios 

Because of its excellent particle identification 
capabilities, MIPP upgrade will measure the ra- 
tio of charged kaons to pions as a function of 
Pl,Pt of the final state particle. This measure- 
ment is of importance to both the NuMI tar- 
get measurements and the atmospheric neutrino 
measurements, since the charged K's produce v' e s 
which are a background to the oscillation search 

Vyi -> Ve- 

2. MIPP Upgrade Proposal 

In October 2006, MIPP proposed to upgrade its 
detector by increasing the data acquisition speed 
of its TPC from 20 Hz to 3000 Hz using the AL- 
TRO/PASA chips developed for the ALICE col- 
laboration [17J. With this large factor in data 
acquisition speed, it becomes possible to acquire 
5 million events in a single day of running. The 
electronics of the rest of the detector will also be 
upgraded to run at this speed. We propose to use 
a hemisphere of the plastic ball detector to mea- 
sure the target recoil particles. This enables us 
to measure particle production on 30 of the most 
common nuclei found in particle physics detec- 
tors and improve the hadronic shower simulator 
programs. It will also permit us to study non- 
perturbative QCD in unprecedented detail. The 
baryon resonance spectrum can be investigated 
up to 3 GeV/c 2 , using pion and kaon beams in 
the momentum range 1 GeV/c-5 GeV/c. The de- 
cision on the upgrade proposal has been deferred 
till analysis results from the present run are avail- 
able and the collaboration has been strengthened 
further. 

2.1. Upgrade Run Plan 

We propose to conduct the upgrade running 
in three phases. In the first phase, we will ac- 
quire data on the neutrino targets, liquid nitrogen 
and take particle production measurements on 12 
thin nuclear targets (5 million events per nucleus) 
which will improve the quality of hadronic shower 
simulation programs. The number of events ac- 
quired in each target is indicated in table [U 

During phase 2, we plan to complete the re- 
maining 18 nuclei of the A-List, as detailed in the 



Target 


Events 


Running Time 




(Millions) 


(Days) 


NuMI target 1 


10 


2 


NuMI target 2 


10 


2 


Liquid Hydrogen 


20 


4 


Liquid Nitrogen 


10 


2 


12 Nuclei 






D 2 ,Be,C,Al,Si,Hg, 






Fe,Ni,Cu,Zn,W,Pb 


60 


12 


Total Events 


110 


22 



Table 1 

Phase 1 Run Plan. 



Target 


Events 


Running Time 




(Millions) 


(Days) 


18 Nuclei 






Li,B,0 2 ,Mg,P, 






S,Ar,K,Ca,Ni,Nb,Ag, 






Sn,Pt,Au,Pb,Bi,U 


90 


18 


10 Nuclei B-list 






Na,Ti,V, Cr,Mn,Mo, 






I, Cd, Cs, Ba 


50 


10 


Total Events 


140 


28 



Table 2 

Phase 2 Run Plan. 



section on hadronic shower simulation and then 
proceed with the B-list if there is need. The sec- 
ond phase of running is detailed in table El 

During phase 3, we plan to go into the tagged 
neutral beam mode, where we run the liquid 
hydrogen target and allow the ILC calorime- 
try to run simultaneously in place of the MIPP 
calorimeter to study the response of the ILC 
calorimeters to tagged neutral beams. 

3. Conclusions 

MIPP has acquired high quality particle pro- 
duction data which it is busy analyzing. The 
MIPP upgrade will improve quality and statis- 
tics data by an order of magnitude by speeding 
up the upgrade. Collaborators are welcome. 
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